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ABSTRACT 

Massive stars at redshifts z .> 6 are predicted to have played a pivotal role in cosmological reion- 
ization as luminous sources of ultra-violet (UV) photons. However, the remnants of these massive 
stars could be equally important as X-ray luminous (Lx ~ 10 38 erg s _1 ) high- mass X-ray binaries 
(HMXBs). Because the absorption cross section of neutral hydrogen decreases sharply with photon 
energy (a oc E~ 3 ), X-rays can escape more freely than UV photons from the star-forming regions in 
which they are produced, allowing HMXBs to make a potentially significant contribution to the ioniz- 
ing X-ray background during reionization. In this paper, we explore the ionizing power of HMXBs at 
redshifts z ;> 6 using a Monte Carlo model for a coeval stellar population of main sequence stars and 
HMXBs. Using the archetypal Galactic HMXB Cygnus X-l as our template, we propose a composite 
HMXB spectral energy distribution consisting of black-body and power-law components, whose con- 
tributions depend on the accretion state of the system. We determine the time-dependent ionizing 
power of a combined population of UV-luminous stars and X-ray luminous HMXBs, and deduce fit- 
ting formulae for the boost in the population's ionizing power arising from HMXBs; these fits allow 
for simple implementation of HMXB feedback in numerical simulations. Based on this analysis, we 
estimate the contribution of high rcdshift HMXBs to the present-day soft X-ray background, and 
we show that it is a factor of ~ 100 — 1000 smaller than the observed limit. Finally, we discuss the 
implications of our results for the role of HMXBs in reionization and in high redshift galaxy formation. 
Subject headings: galaxies: formation - X-rays: binaries - cosmology:theory 



1. INTRODUCTION 

There is strong and compelling observational evidence 
that the Universe underwent an "Epoch of Reioniza- 
tion" within__tb£ first_~ l_biJliojiyeaxsafter_rtieJBiK 
fe.g.lOuchi et al.ll2010t IMesingenl20Tol IShull et al.l[20ll 
iMcGreer et alj 120111 ). During this period, the cosmic 
abundance of neutral hydrogen was "re-ionized" by a 
background of ionizing ultra-violet (UV) and X-ray radi- 
ation pr oduced by the first gene r ation of stars and galax - 
ies (e.g. IBarkana fc Loebl 120071 IRobertson et all 120101 ) . 
The precise nature of these sources remains an outstand- 
ing problem, but it can be argued reasonably that mas- 
sive stars (M* ;> 8M Q ) must have been important (e.g. 
IWise fc Abell2lM IWIs3l20T2h . 

Massive stars have fleeting main sequence (MS) 
lives (~10 Myrs), but as extremely luminous sources 
of hydrogen-io nizing UV photons (see, for example, 
iSchaererl 120031 ) they are expecte d to play a crucial role 
m rciomzmg the Uni verse (e.g. iWvithe fc Loebl 120031: 
ISokasian et all l200i IWise fc Abe J I2008D . However, 
their predicted contribution to reionization depends 
on the ease with which UV photons can escape into 
the intcrgalactic medium (IGM), which depends on 
vario us factors, principa lly the dumpiness of the IGM 
(e.g. iPawlik et al.1 l2009t) . However, testing such pre- 
dictions observationally is inherently challenging and a 
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clear consensus as to what the sources of reionization 
are is yet to emerge. For example, observations of 
star-forming galaxies at redshifts 7 <J z <J 10 hint that 
implied star formation rates are insufficient to produce 
enough massive stars to reionize the Universe by z ~ 6, 
without including a population of galaxies below the 
detection limits and relaxing assumptions about the 
escape fraction of ioni zing photons and ex t rapolated star 
form a tion rates (e.g. IBunker et al.l 120101: iMcLure et al.l 
120101 : IFinkelstein et all 120101 : iGonzalez et al.l 120101 ). Ln 
contrast, observations of high redshift galaxies also 
suggest that the star formatio n rate at z ~ 9 could 
be as high as that at 2=0 (cf. llshida et all 1201 If ) and 
that the observable population of galaxies between 
6 <J z <J 8 can sustain a fully reionized IGM at z~6 if 
the a verage escape fraction of ionizing p hotons is ~30% 
(cf. IFinkelstein et all l2012t see also iBouwens et al 
I2012D . A similar conclusion is drawn by iKistler et al 
(|2009l) . who argue that a sufficient number of mas- 
sive stars could have formed to reionize the Universe 
by using high-redshift gamma ray bursts as a proxy 
for the integrated star formation rate between 6 <J z <J 8. 

It is, however, worth noting that massive stars can 
continue to ionize during the ir post-MS lives, predom- 
inantly as sources of X-rays. lOhl (|200ll) noted that X- 
ray emission from supernovae in star forming regions in 
the high-redshift Universe should be large and compara- 
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ble energetically to UV emis sion. The role of supernovae 
has been explored further bv lJohnson fc Khochfarl ([201 ID 
who examined how strong shocking of the ISM associated 
with supernovae leads to the production of ionizing pho- 
tons with harder spectra and larger escape fractions than 
UV photons. They estimate that such X-rays can boost 

- briefly - the ioniz i ng pow er of a massive star by ~ 10%. 
iGlover fc Brand! ([20031 ) considered not only X-rays 

from supernovae but also from X-ray binaries in their 
examination of the contribution of star formation to 
the build up of the high re dshift X- ray background. 
As Uustham fc Schawinskil ([2012D have demonstrated 
in their recent study, X-ray binaries are likely to be 
imp ortant sources of fe edback across cosmic time (see 
also iFragos et al.l l2012t ) . In partic ular, we expect this 
to be the case at early times; in IPower et al.l (|2009f ) 
(hereafter P09) we explored how high mass X-ray 
binaries (HMXBsfl in primordial globular clusters at 
high redshifts could boost the UV ionizing power of the 
cluster. We found that harder ionizing spectra combined 
with enhanced escape fractions for X-rays implied that 
HMXBs could be just as efficient at ionizing the IGM as 
th eir MS progen i tors. Simil ar ideas have been explore d 
in lMirabel et al.l ([20111 ) and lWheeler fc Johnson] (|2011| ). 
in wh ich the primary i s a stellar mass black hole, as well 
as by [McQuinn] ([20121 ) in his census of potential sources 
of reionization. 

The analysis in P09 focused on the ionizing power of 
HMXBs in globular clusters because (i) the inferred ages 
of metal p oor globular clusters imp ly that they formed at 
z £ 6 (cf. iBrodie fc Straderl [20061) ; (ii) the relationship 
between the initial mass functions (IMF) of stars and 
the dynamical evolution of c lusters is well understood 
(cf. IVesperini fc Heggid Il997l ) : and (hi) the escape frac- 
tion for UV photons was likely to be large, assuming that 
globular clusters fo llowed similar orbits then to the ones 
they follow now (cf. Ricotti 2002, who explored the ioniz- 
ing power of massive stars in globular clusters). This last 
point is important because the large UV escape fraction 
allowed P09 to carry out a straightforward comparison 
of a globular cluster's UV and X-ray ionizing power. 

However, as noted in P09, HMXBs are likely to be 
a generic by-product o f high mass star formation (cf. 
iHelfand fc Moranl I2001D . and should form as well in 
gas-rich galaxy discs as in young globular clusters. This 
is consistent with observations of star forming galaxies 
at z <^ 1 that show that luminous compact X-ray sources 

- with properties mirroring those of HMXBs in our 
Galaxy - are good indic ators of recent star formation 
activ ity (see Figure 10 of iMineo et al.ll2012al and Figure 
2 of IMineo et al.1 l2012b| ). In this paper, we extend 
the analysis presented in P09 to study HMXBs as 
generic sources of ionizing radiation. Here our focus 
is on refining our treatment of their spectral energy 
distribution. Rather than assuming a simple power-law 
form for the HMXB spectrum (as we did in P09), we use 
the spectrum of Cygnus X-l, the archetypa l Galactic 
HMXB (e.g. iRemillard fc McClintockl 120061 ). as our 
template. We assess how this impacts on the ionizing 
power of a coeval stellar population over time and we 

* We assumed an initial primary mass of M* > 8Mq . 



quantify how the ionizing luminosity of the population 
is boosted by the presence of HMXBs. Finally, we 
estimate the possible contribution of HMXBs to the 
present-day soft X-ray background, assuming both our 
new template spectrum and power-law spectra of the 
kind that have been used in previous studies (e.g. P09). 

The structure of the remainder of this paper is as fol- 
lows. In §[2] we describe our time-dependent Monte Carlo 
model for the spectral energy distribution of a coeval 
population of stars in which HMXBs are forming. In § [3] 
we present results for the time evolution of the ionizing 
power and spectral energy distribution of the population, 
and we quantify how HMXBs enhance the ionizing power 
of the stellar population. In § [4] we show that the X-ray 
luminosity produced in our model does not violate ob- 
served limits on the soft X-ray background. Finally, we 
summarize our results in § [5] and comment on their im- 
plications for cosmological reionization and high redshift 
galaxy formation. 

2. METHODS 

2.1. Modelling HMXBs in a Single Stellar Population 

As in P09, we set up a Monte Carlo model of a stellar 
population, assumed to form in a single instantaneous 
burst, and follow the evolution of the massive stars over 
the first 250 million years, through their MS lives and 
into the HMXB phase. The main features of our model 
can be summarized as follows; 

THE MASSIVE STAR POPULATION: We assume the 
IMF of lKroupal pOOll ) as our fiducial case, with stellar 
masses spanning the range 0.0 1 < M^/Mp, < 100, but 
we als o verify our resu lts for the lSalpeterl (|1955|) and top- 
heavy [Ch^riei (|2001[ ) IMFs. All stars with M* > 8M Q 
are assumed to form in binaries - these are the progen- 
itor population from which the HMXBs are drawn. Ini - 
tial binary parameters are assigned following iDravl ()2006l ) 

- that is, companion masses are drawn from a uniform 
distribution between 0.01 < M*/M < 100 and orbital 
periods are distributed uniformly in logarithm between 1 
and 10 4 days. Massive star lifetimes are e stimated using 
the re sults of iMarigo et all (I2001L I2003D . iSchaerer et "all 
(|1993f ) and lMevnet fc Maederl (|2000l ) for metallicities of 
Z = 0, 0.008 and 0.02 (i.e. solar metallicity) respectively; 
we explore the Z = case in our results section. 

THE HMXB POPULATION: We assume that HMXBs 
form from binaries in which the initial MS mass of the 
primary exceeds M* ~ 8Mq, the threshold for neu - 
tron star formation (cf. Figure 1 of lHeger et all[2003al) . 
and the donor (i.e. secondary) mass lies in the range 
> 3Mg0. Once the primary goes super nova, we es- 
timate the remnant mass using Figure 3 of lHeger et al.1 
(|2003al) . The binary will be disrupted if it loses more 
than half its mass in the supernova - which, for our fidu- 
cial Kroupa IMF, implies that approximately 30% of bi- 

t This donor mass is lower than t he usual definiti on of HMXBs 

- donor OB stars (cf Table 1 of Hj bbiano 2006) with typical 
masses > IOMq (e.g. Justham & Schawinski 2012). However, it 
is reasonable to include these intermediate mass X-ray binaries 
because they are sufficiently luminous to contribute to the X-ray 
ionizing power of the stellar population and their MS lifetimes are 
of order ~ 10 s yrs. 
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naries survive, thereby setting an upper limit to the total 
number of HMXBs that can potentially form. Follow- 
ing P09, we draw a survival fraction / sur of this ~ 30% 
at random and consider them as HMXBs; / sur captures 
the various uncertainties that prevent massive binaries 
from evolving into HMXBs. We assume that HMXBs 
are active until the companion star evolves off the main 
sequence and goes supernova. 

We have noted already that the binaries that survive 
are more likely to host black holes than neutron stars, 
especially in low metallicity systems where the forma- 
tion rate of HMXBs in which the primary is a black hole 
could be a factor of ~ 10 higher than at Solar metallic - 
ity (cf. lLinden et alj|2010t [Justham fc Schawinskill2012D . 
There is also observational evidence that the black hole 
mass is likely to be larger in lower metallicity sys- 
tems (cf. iCrowther et al.l l2010f ). This should, m prm- 



Table 1 

Adopted Mass-Luminosity Relation. 



ciple, shape our HMXB luminosity function (cf. iDravl 
120061) . However, we make the simplifying assumption 
to draw HMXB luminosities from a Weibull distribu- 
tion with a peak luminosity of Lx ~ 10 38 erg s" 1 
but capped such that they do not accrete at super- 
Eddington rates; this sets an upper limit of approxi- 
mately L x ^ 1.26 x 10 38 (M/M Q ) erg s" 1 on the luminos- 
ity of an HMXB with primary mass M. This approach 
gives a distribution that is consistent with the luminosi- 
ties of compact X-ray sources in nearby galaxies whose 
X-ray binary populations are dom inated by HMXBs (cf. 
Figure 1 of iGilfanov et "all 120041 ) : see P09 for further 
discussion of this point. 

2.2. Spectral Energy Distribution & Time Dependence 

We split the spectral energy distribution of the stellar 
population into two components: 

MAIN SEQUENCE (MS) STARS: Each star is assumed to 
radiate as a black body with an effective temperature of 



sb) 



1/4 



(1) 



where R*/Rq = (M*/M0) - 8 is stellar radius and ctsb is 
the Stcfan-Boltzmann constant. Stellar luminosity L. t is 
assumed to follow a mass-luminosity relation of the form 



M* 







(2) 



where a, which governs the amplitude of the relation, 
and /?, which go verns its slop e, are observed to depend on 
stellar mass (e.g. lHenrvll2004l ): their values and the stellar 
mass range in which they are applicable are summarized 
in Table III 

Empirically there is a wealth of evidence that the 
rate at which stellar luminosity varies with mass (J3 = 
d log L/d log M) d e crease s with increasing mass M*. For 
example, IMalkovi (|2007f ) find j3 ~ 4.1 at M ~ 1 M Q 
to ~ 3.2 at M ~ 20 M Q (see their Table 6), while 
IVitrichenko et al.l (|2007l) report that /3 ~ 2.76 over the 
mass range 20 < A//M© <J 50. Our values of a and (3 
provide a good approxim ation to the functional form pre- 
sented in IMalkovi (|2007l ) for stellar masses M* £ 50M o . 
Above 50M Q we must extrapolate because observation- 
ally inferred data are few; we assume that the relation 
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is slightly shallower than it is below 50Mq. This is an 
uncertainty, but it has negligible effect on our results be- 
cause relatively few stars are formed in this mass range 
and their lifetimes arc short. We find that the mass- 
luminosity relation can be well approximated by a 2 nd - 
ordcr polynomial of the form, 

logL»/L Q ^co + CilogM^/Mo+caQogM./Mo) 2 . (3) 

Here the coefficients have the values Co=-0. 04172, 
ci=4.4954 and C2=-0.6041, and log indicates logarithm 
base 10. 

HIGH MASS X-RAY BINARIES: We model HMXB spec- 
tra using the archetypal Galactic HMXB Cygnus X- 
1 as our template. Cygnus X-l is the brightest 
HMXB in the Galaxy and it has been studied in 
exquisite detail (see, for exam ple, the recent review by 
iRemillard fc McClintockl 120061). It consists of a black 
hole of mass 8.7±0.8M g (|Shaposhnikov fc Titarchukl 
120071) and a super-giant companion (HDE 226868). Its 
spectrum is observed to fluctuate between distinct low- 
hard an d high-soft stat e s, exa mples of which can be 
found in lGierlihski et al.l (|1999l ). 

• The low-hard state is characterized by a luminosity 
of -L2-iokcV ~ 3xl0 36 erg s _1 and a hard power- 
law spectru m (xE~ r with spectral ind ex of T ~ 
1.6-1.8 (e.g IGierlihski etlllll997LIT99l . 

• The high-soft state is characterized by luminosity 
roughly one order of magnitude higher and a strong 
black body component with kT ~ 0.5 keV with soft 
power-law tail with T ~ 2.5 fe.g. lDolan et al.lll977l: 
lOgawara et al.lll982] ). 

This suggests that we should adopt distinct spectral 
shapes for low-hard and high-soft states. However, we 
understand neither what sets the duration of these two 
states - for example, Cygnus X-l is observed to be pre- 
dominantly in its low hard state, but an alogous systems 
such as LMC X-3 (|Val-Baker et al.l2007l ) appear to spend 
more time in their high soft state - nor how the na- 
ture and duration of these states depend on factors such 
as metallicity. For this reason, we make the simplify- 
ing assumption that HMXB spectra do not vary in time 
and instead introduce a threshold in X-ray luminosity of 
10 37 erg s _1 , above (below) which a source is in a soft 
(hard) state. In particular, we model 

• the low-hard state as a power-law of slope -0.8 be- 
tween 2 — 10 keV, the energy range in which this 
spectrum profile is observed, and 

• the high-soft state by a composite black body and 
power-law spectrum - the power-law has a slope 
-1.1 and is normalized such that the ratio between 
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Figure 1. Spectra of the total energy output of the fiducial stel- 
lar population (i.e. N=10 6 , Kroupa IMF, / sur = 1) after 5, 50, 
100, and 150 Myrs (solid, long-, medium- and short-dashed curves 
respectively). The spectral shape derives from the stellar UV black 
body and the HMXB softer black body + harder power-law com- 
ponent, where we use the spectrum of Cygnus X-l for our HMXB 
template. 

the luminosities of the black body and power-law 
components matches that of Cygnus X-l, while the 
black body temperature is calculated by assuming 
(Lx/Lcy g x-i) = (T/Tcygx-i) 4 - 

3. RESULTS 

In the following subsections, we show how the ioniz- 
ing power of a stellar population formed in an instan- 
taneous burst evolves over the first 150 Myrs after for- 
mation (§ 13. ip . We compare and contrast results in the 
presence and absence of HMXBs and, where appropriate, 
we comment on the sensitivity of our results to the as- 
sumed value of the survival fraction (/ SU r) and our choice 
of IMF. 

3.1. Ionizing Power over Time 

In Figure [T] we show the spectral energy distribution of 
the fiducial stellar population (i.e. Kroupa IMF, / sur = 1 
and N=10 6 ) as a function of time - a composite of MS 
stars, which we model individually as black bodies, and 
HMXBs, which we model as a combination of black bod- 
ies and power-law components. The amplitude of the 
black body corresponding to MS stars decreases with 
time, while its peak shifts to lower energies; this reflects 
the evolution of the most massive stars in the stellar pop- 
ulation - which dominate the UV-luminosity - off the 
MS. Over the same period the amplitudes of both the 
black body and power-law components of the HMXBs 
decrease with time. 

The same qualitative trends can be seen if we vary the 
IMF from Kroupa to Chabrier or Salpeter - there is a 
systematic increase (decrease) in the amplitude for the 
Chabrier (Salpeter) IMF, which reflects an increase (de- 
crease) in the proportion of massive stars that form. If 
we vary the size of population between N=10 4 to iV=10 6 
the amplitude varies linearly with the size of population. 
Varying the survival fraction / sur between 0.01 and 1 
suppresses the amplitude of the HMXB black body and 
power- law contribution while leaving the black body 
contribution from the MS unchanged. 

In Figure [2] we estimate how much of this energy (in- 
stantaneous in the upper panel, cumulative in the lower 
panel) is available to ionize neutral atomic hydrogen (HI) 
in the IGM as a function of time. We evaluate this by in- 
tegrating the spectra plotted in Figure [T] between a lower 
limit of 13.6 eV, the minimum photon energy required to 
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Figure 2. Time evolution of the instantaneous and cumulative 
ionizing power of a stellar population of N = 10 6 stars and a 
Kroupa IMF (upper and lower panels respectively). Solid curves 
correspond to the contribution from MS stars only; curves of dif- 
ferent line types correspond to the combined contribution of MS 
stars and HMXBs for five different values of the survival fraction 

/sur with 0.01 < /sur < 1. 

ionize atomic hydrogen, and an upper limit of £ii m , the 
energy above which atomic hydrogen becomes transpar- 
ent to hard X-ray^3; for the redshift range we consider 
(z £ 6), E hm £ 1 keV. Provided / sur > 0.1, we find 
that HMXBs dominate the ionizing power of the stellar 
population after ~ 20 Myrs. 

Interestingly we note that the X-ray ionizing power 
shows a pronounced bump between 20 <; t <J 40 Myrs. 
This is an imprint of our definition of an HMXB as a 
binary system in which the mass of the primary must ex- 
ceed ~ 8M Q at the end of its MS lifetime - the threshold 
for neutron star formatio n. Because the MS lifetim e of 
an 8M is - 25 Myrs (cf. iMarigo et aill200l l200l and 
because statistically we expect more HMXBs to be drawn 
from stars close to this 8M Q limit, we expect an excess 
in HMXBs to form around this time and to produce a 
bump in the ionizing power. After this time, HMXBs 
can no longer form and the ionizing power comes from a 
population that is in terminal decline. 

We quantify how HMXBs enhance the ionizing power 
of the stellar population by measuring the boost factor 
/hmxb = Lhmxb/Lms, which compares the ionizing lu- 
minosity of the MS stars only (Lms) an d of the com- 
bined MS stars and HMXBs (Lhmxb)- In Figure [3] we 
show how /hmxb varies with time for a stellar population 
of 10 6 stars and / SU r=l assuming Kroupa, Salpeter and 
Chabrier IMFs respectively (solid curves, left to right 
panels). Because the ionizing power of the MS compo- 
nent declines sharply after ~ 10 — 20 Myrs as the most 
massive stars end their lives, /hmxb grows steadily with 

' As in P09, we estimate Eii m by requiring that <r(i?ii m ) = 
c~ 1 H(z) njj(z) -1 where cr(E) is the ionization cross section of neu- 
tral hydrogen, H(z) is the Hubble parameter at redshift z, nn(z) 
is the mean number density of hydrogen and c is the speed of light. 
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Figure 3. Time evolution of the boost factor /hmxb = 
-^hmxb/^ms f° r a stellar population of N = 10 6 stars, a sur- 
vival fraction / sur = 1 and a Kroupa, Salpeter and Chabrier IMFs 
(left, middle and right panels respectively). Dotted-dashed curves 
correspond to polynomial fits of the form given by Equation[4]with 
parameters given in Table [2] 



Table 2 

Boost Factor /hmxb Parameters. 



IMF 



Kroupa 7.85 
Salpeter 6.53 
Chabrier 0.852 



6.73 
6.82 
8.23 



54.2 
66.2 
46.8 



-85.4 45.0 
-104 55.1 
-69.2 34.3 



-7.61 
-9.40 
-5.42 



time, peaking at /hmxb ~ 1600 after t ~ 180 Myrs be- 
fore dropping off abruptly. The dashed curves in Figure[3] 
show fits to the numerical results that agree to better 
than 10%. We fit the results with 2 formulae, below and 
above the bump at 30 Myr. These are given by 



log(/] 



HMXB ) = 



10 -5 p e .3iogt 
a + blogt + c(logt) 2 



t < 30 



(4) 



-d(logtyt > 30 



where t is population age in Myr, and are valid up to ap- 
proximately 180 Myr, above which HMXB contribution 
falls rapidly to (and /hmxb drops to 1). The parame- 
ters p, q and a to d are tabulated in Table [3] 

4. COSMIC SOFT X-RAY BACKGROUND 
The Universe at z > 6 is effectively transparent 
to photons with energies ;> 1 keV. This means that 
hard photons remain unabsorbed by neutral hydrogen 
and so a hard X-ray background was built up in the 
early Universe that, through redshifting, contributes 
to the present day soft X-ray background (SXRB), 
the mean specific background intensity of soft X- 
rays. We wish to determine the contribution made 
to the SXRB by the HMXB population present at 
z > 6 and whether or not their contribution violates 
observed limit s. Equivalent calcu lat ions have been 
perfor med bv iDiikstra et ahl (120041) , iSalvaterra et al.1 
pfM l and iZaroubi et all (|2007l ) for mini-quasars in 
th e early Universe. Ho wever, we adopt the approach 
of IDiikstra et al.l ([2012D , who looked at the contribu- 
tion of star-forming galaxies to both the soft and hard 
X-ray backgrounds from high redshift to the present day. 

In what follows, we assume values of Hq = 
100ft. km s _1 Mpc -1 for the Hubble parameter with 
ft=0.71, f^M = 0.266 for the matter density param- 
eter, and S1a = 0.734 for the dark energy den- 
sity parameter (cf . iKomatsu et al.l 120 111 ) , and we take 
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Figure 4. Comparison of the composite (i.e. black-body plus 
power-law) HMXB spectral energy distribution (solid curve) and 
possible power-law alternatives (with indices a=0.5,1.5,2.5, corre- 
sponding to dotted, dashed and dottcd-dashed curves respectively). 
The composite HMXB spectrum is as measured at i=50 Myrs in 
Figurc[T] while the power-law spectra are normalized such that the 
total energy emitted between 1 and 2 keV is the same for both the 
composite and power-law spectra. 



the 



SXRB the 

2 



value of 



as our observed limit 

~ 3.4 ± 1.4 x 10~ 1 3 ergs ~ 1 cm~ 1 deg~ z obtained by 
iHickox fe Markevitchl ([2007D for the excess flux between 
1 — 2 keV, which corresponds to hard X-rays at z ;> 6. We 
estimate the contribution made by high redshift HMXBs 
to the observed SXRB by evaluating 



SXRB 



An c 



p*£x(z) 
(l + z) 2 S(zY 



-dz, 



(5) 



where AQ ~ 3 x 10 4 srdeg 2 ; c is the speed of light; 
z max is the earliest (unknown) redshift at which HMXB 

formation proceeds; £(z) = a/JIm^ + z) 3 + £1a; P* is 
the comoving star formation rate density in units of 
M Q yr -1 Mpc~ 3 ; and Cx(z) is the "K-corrected" X-ray 
luminosity per unit star formation rate in the observed 
energy range X = E\ — E 2 ■ 

For we explore four possibilities, spanning the range 
of potential star formation histories. We take (i) a fixed 
value of p*(z = 6) c± 0.002 M^ yr" 1 Mpc -1 , as mea- 



sured at z=6 bv iBunker et al.l (|2010D ; (ii) a fixed value 
of p*(z — 6) ~ 0.17 Mf^ yr" 1 Mpc" 1 , as measured 
at z=6 by iTanvir et aD (|2012l ): (iii) a fixed value of 



p*(z = 6) ~ 0.02 M 



Mpc" 



estimated at z=6 



using the functional form of iHopkins fc Beacoml (|2006| ). 



P*{z) 



(u + vz)h 
l + (z/w) x ' 



(O 



with parameters m=0. 017, f=0.13, u>=3.3 and £=5.3 
from iCole et~aTI (|200l ; and (iv) a rcdshift-dcpcndcnt 
value of p* , estimated using Equation [6] 

For £x(z) we follow lDiikstra et al.l (|2012f ) and recast it 
as Cx(z) — cxK(z). Here we assume the range of values 
for 2.6 < cx / lO^erg s-HM^yr- 1 )- 1 < 3.7 measured by 
iMineo et al.l (|2012bl ) for compact resolved X-ray sources 
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in galaxies (lower limit) and unresolved galaxies in the 
Chandra Deep Field North and ultra-luminous infra-red 
galaxies (upper limit); and we calculate 



K(z) 



E 2 (l+z) 
£1(1+2) 



EF(E)dE 



(7) 



where F(E) represents the spectrum of a single HMXB, 
Ei = l keV and i?2=2 keV. Examples of the spectra we 
consider are shown in Figure [4] - the solid curve corre- 
sponds to the cumulative spectrum derived from HMXBs 
in the fiducial stellar population described in §3, as 
measured at i=50 Myrqa, while the dashed, dotted and 
dottcd-dashed curves correspond to power-law spectra 
(oc E~ a ) that have been normalized such that the to- 
tal energy emitted between 1 and 2 keV is the same for 
both composite and power-law spectra. The units are 
arbitrary - it is the shape, not the amplitude, that is 
important when evaluating Equation [7] - and it is worth 
noting how poorly simple power-law spectra approximate 
the shape of the more realistic composite spectrum. 

In all cases we fix our lower bound in Equation [5] at 
z — 6 and allow our upper bound z max to vary between 
z = 20 and z=50. This corresponds to a range in age 
of the Universe of between 48 and ~ 180 Myrs old 
(z=50 and 20 respectively), which is likely to bracket 
the redshifts when the first generation of dark matter 
halos were sufficiently massive to support cooling and 
subsequent star formation formed (cf. Figure 1 of 
Glover 2005). However, we note that our results are in- 
sensitive to z max for the range of values that we consider. 

Evaluating Equation ([5]), we find that the contribu- 
tion of HMXBs to the soft X-ray background is of order 
~ 5 x 10~ 16 erg s _1 cm~ 2 deg~ 2 , depending on what 
precisely is assumed for the star formation rate density 
and spectrum - roughly a factor of ~ 1000 smaller than 
the observed limit of ~ 3.4 x 10~ 13 erg s" 1 cm -2 deg -2 
measured by iHickox fc Markevitchl (|2007fl . In general, 
the higher the star formation rate density, the larger the 
contribution that is possible, because of our assumption 
that HMXB formation tracks (massive) star formation. 
If power-law spectra are adopted, we find that SXRB 
~ 2.15 x 10~ 16 erg s _1 cm -2 deg~ 2 for a=0 compared to 
~ 11.4 x 10~ 16 erg s _1 cm~ 2 dcg~ 2 for a=2; that is, the 
softer the spectrum, the larger the contribution HMXBs 
can make. Nevertheless, this contribution is still a factor 
of ~ 100 smaller than the observed limit. 

5. SUMMARY 

There are compelling astrophysical reasons to expect 
that HMXBs, which are observed to be a natural by- 
produ c t of massive star fo rmation (e.g. iMineo et al.l 
l2012at iDijkstra et all l20ll . could be an important 
sourc e of feedback in galaxy forma tion over cosmic time 
(e.g. Uustham fc Schawinskil l2012f ). Because the cross 
section for photons to be absorbed by neutral hydro- 
gen decreases strongly with energy as E~ 3 , X-rays from 
HMXBs can effectively diffuse through the IGM and de- 
posit their energy at relatively large distances from the 

§ We treat the unnormalized cumulative spectrum as indicative 
of the shape of a typical HMXB. 



stellar population, in an act of non-local feedback (a 
similar effect is expected fo r mini-quasars; see, for ex- 
ample, [Zaroubr^Silk 2005). This contrasts with lower 
energy UV photons, whose correspondingly larger ab- 
sorption cross section means that they ionize the IGM in 
the vicinity of the stellar population, in an act of local 
feedback. 

In t his paper, we build on the work of iPower et al.l 
(|2009t ) in which we examined the X-ray luminosity and 
effective ionizing poweid of a coeval population of stars 
and HMXBs at z ;> 6. 

• We have used the archetypal Galactic HMXB, 
Cygnus X-l, as our template for a more realistic 
HMXB spectrum. The form of this spectrum de- 
pends on the accretion state of the system and is 
a composite of black-body and power-law compo- 
nents. 

• Using this new composite spectrum, we have quan- 
tified how HMXBs enhance the ionizing power of 
the stellar population by defining the boost factor 
/hmxb = Lhmxb/Lms (where Lms and Ihmxb 
are the ionizing luminosities of a population of MS 
stars and a combined population of MS stars and 
HMXBs respectively) and characterizing and its 
variation with time. Because the ionizing power of 
the MS component declines sharply after ~ 10-20 
Myrs as the most massive stars end their lives, we 
find that /hmxb grows steadily with time before 
peaking at t ~ 100 Myrs and after dropping off 
abruptly. This demonstrates the relatively long- 
lived nature of HMXBs as ionizing sources when 
compared to young massive stars on the MS. 

• Finally, we have estimated the contribution of 
HMXBs to the soft X-ray background by assum- 
ing that HMXB formation is linked explicitly to 
the global star formation rate density, and we show 
that it does not violate observed limits. Depending 
on what is assumed for the star formation rate den- 
sity, we estimate a contribution of ~ 5 x 10~ 16 erg 
s _1 cm~ 2 dcg~ 2 if HMXB spectra are modelled as 
composite black-bodies and power-laws - roughly 
a factor of ~ 1000 smaller than the observed limit 



of ~ 3.4 x 10 13 erg s 1 cm 2 deg 2 measured by 



IHickox fc Markevitchl (|2007l) . 

We do not compare the HMXB contribution to the 
ionizing X-ray background present during reionization to 
the contribution from other source populations, such as 
mini-quasars and accreting super-massive black holes, in 
this paper. Suc h a comparison has been carried out by 
iMcQuinnl (|2012[ ) who concludes that HMXBs were likely 
to be a minor contributor to reionization. This may be 
the case - although the conclusion depends on the as- 
sumed spectral properties of the various sources consid- 
ered, which in most cases are uncertain - but our work 
implies that HMXBs can have a significant impact on 
the ionization structure and heating rate of the IGM at 

^ We made the simplifying assumption that an energetic X-ray 
photon can ionize multiple hydrogen atoms, and so treated sec- 
ondary electrons that ionize hydrogen atoms as effective photons. 
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z J> 6, extending the volumes of partially ionized hydro- 
gen by factors of ;> 1000. This will have a direct impact 
on the efficiency of galaxy formation by, for example, sup- 
pressing t he collapse of gas onto low-mass dark m atter 
halos fe.g. lMachacek et al.|[2003HTanaka et al.ll2012[ ) and 
modifying the co oling rate of hot gas in galaxy halos (e.g. 
ICantalupoll2010( ). We will demonstrate this explicitly in 
a forthcoming paper (James et al., in preparation), in 
which we will use the more realistic composite spectrum 
in a radiative transfer calculation to assess the relative 
importance of such HMXB feedback in heating the high 
rcdshift IGM and the potential HMXB contribution to 
both cosmological rcionization and galaxy formation in 
general. 
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